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ABSTRACT 


The  work  conducted  at  Louisiana  State  University  on  a  class  of  prob¬ 
lems  arising  in  the  evaluation  of  weapons  and  in  the  analysis  of  the  best 
methods  for  their  use  is  reported.  The  destructiveness  of  a  specified 
weapon  against  a  specified  target  is  measured  by  the  ratio  of  the  number 
of  targets  destroyed  to  the  number  of  weapons  used. 

This  report  describes  the  methods  developed  for  the  analytical  eval¬ 
uation  of  average  kill  probabilities  for  single  weapons  and  proposes  the 
optimization  of  pattern  parameters  for  multi-weapon  attacks  on  discretely 
defined  target  complexes  for  future  investigation. 

An  acceleration  algorithm  for  minimizing  a  convex  objective  function 
subject  to  linear  constraints  is  also  presented  in  detail. 
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SECTION  I 


INTRODUCTION 


This  is  a  report  of  the  work  conducted  at  Louisiana  State  University 
on  a  class  of  problems  arising  in  the  evaluation  of  weapons  and  in  analyz¬ 
ing  tne  best  methods  for  their  use.  The  effectiveness  (destructiveness) 
of  a  weapon  against  a  specified  target  can  be  measured,  at  least  in  part, 
by  the  ratio  of  the  number  of  targets  destroyed  to  the  number  of  weapons 
used.  The  number  of  weapons  required  to  destroy  a  target  depends  on  two 
primary  factors:  the  probability  that  the  target  is  destroyed  if  the 
weapon  hits  it  and  the  probability  that  the  weapon  will  hit  target.  If, 
in  addition,  weapons  are  used  in  groups,  the  results  also  depend  on  the 
firing  pattern  employed.  The  multi-weapon  case  was  not  considered  by  the 
researchers  in  the  work  performed  to  date.  This  report  concentrates  on  the 
methods  developed  for  calculating  the  probabilities  of  destroying  targets 
with  single  weapons. 


SECTION  II 


LETHAL  AREA 

The  simplest  case  of  measuring  the  destructiveness  of  a  weapon  against 
a  target  occurs  when  the  weapon  must  hit  the  target  in  order  to  destroy  it, 
but  always  destroys  the  target  if  it  hits  it.  In  this  case,  the  probability 
that  a  single  weapon  will  destroy  the  target  is  just  the  probability  that 
it  hits  a  certain  area,  called  the  Lethal  Area.  The  size  of  the  lethal 
area  is,  consequently,  a  measure  of  the  destructiveness  of  the  weapon 
against  the  target.  If  the  weapon  has  a  proximity  fuse,  it  is  no  longer 
necessary  to  hit  the  target  to  destroy  it,  for  the  proximity  fuse  will 
convert  a  near  miss  into  a  destructive  explosion.  If  the  charge  is  suffi¬ 
cient  so  that  an  explosion  within  the  radius  of  action  R  of  the  proximity 
fuse  will  destroy  the  target,  then  the  lethal  area  is  increased  to  the  area 
included  within  a  curve  surrounding  the  cross  section  of  the  target  at  a 
distance  R  from  its  boundary.  In  the  three-dimensional  case,  one  would 
not  have  a  lethal  area  defined  in  this  manner,  but  a  lethal  volume  instead. 

If  a  single  hit  is  not  enough  to  destroy  a  target,  the  probabilities 
P^  of  destroying  the  target  if  it  is  hit  i  =  1,  2,  ...,  N  times  with  prob¬ 
abilities  must  be  determined.  In  this  case,  the  probability  of  de¬ 
stroying  the  target  is 


P 


k 


N 

=  I 
i=l 


P.P. 
1  i 


CII-1) 


The  probabilities  are  sometimes  called  damage  coefficients.  If  all  the 

damage  coefficients  are  equal  to  each  other,  and  the  probability  of  a  hit 
is  F  for  all  firings,  then,  the  probability  of  destroying  the  target  with 
N  weapons  is,  simply, 

Pk  =  1  -  (1  -  PF)N  (H_2, 

In  other  words,  the  probability  of  destroying  the  target  with  any  one  hit 
is  always  the  same,  regardless  of  how  many  previous  hits  have  occurred. 
This  may  be  interpreted  as  meaning  that  only  when  a  vital  spot  is  hit  will 
the  target  be  destroyed  and  that  hits  elsewhere  on  the  target  will  only 
damage,  but  not  destroy,  the  target.  This  vital  spot  hypothesis  serves 
to  reduce  the  number  of  unknown  quantities  P^  and  F^  to  two,  P  and  F,  and 

has  been  found  to  give  satisfactory  results  in  many  cases,  such  as  AA  hits 
on  aircraft.  When  the  vital  spot  theory  can  be  applied,  the  lethal  area 
of  a  target  is  defined  as  the  product  of  the  effective  area  of  the  target 
and  the  probability  P  that  a  hit  on  this  area  will  destroy  the  target.  For 
example,  the  probability  of  sinking  a  merchant  vessel  with  a  torpedo  hit 
is  about  1/3.  Consequently,  P  =  1/3  and  the  lethal  area  is  1/3  of  the 
ship's  length. 
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In  real  life  situations  it  is  necessary  to  consider  the  variation  of 
the  probability  of  destroying  the  target  as  a  function  of  the  coordinates 
of  the  point  at  which  the  hit  is  made.  In  such  instances,  the  damage 
probability  becomes  a  damage  function  p^Cx.y)  where  x  and  y  are  coordi¬ 
nates  centered  on  the  target  in  a  plane  normal  to  the  weapon's  trajectory. 
In  this  case,  the  lethal  area  may  be  computed  as 


L  =  JJ  P  (x,y)  dxdy . 


(II-3) 


where  the  integration  is  over  all  the  area  for  which  p^  (x,  y)  >  0. 
Formulas  for  one-dimensional  targets  are  obviously  derived  from  (II-3). 

Random  Bombardment 

The  lethal  area  serves  to  evaluate  the  destruction  level  attained  by 
weapons  delivered  at  random  over  an  area  A.  For  any  given  target  within 
the  area  A,  the  probability  that  a  given  weapon  will  hit  the  element  of 
a-ea  dxdy  is  simply  dxdy/A.  The  probability  that  this  weapon  will  destroy 
:he  target  is  therefore 


(1) 

k 


JJ  Pt  (x>y) 


L 

A’ 


(II-4) 


If  N  weapons  are  fired,  the  probability  that  a  given  target  is  destroyed 
is 


(II-5) 


However , 


1/x 

e  =  lim  (1  +  x) 
x-0 


(H-6) 


Therefore,  as  the  number  of  weapons  increases,  i.e., 

N  -  %  (II-7) 


one  can  show  that. 


(N)  ^  -NL/A 
rk 


( II-8) 


This  represents  the  expected  fraction  of  all  targets  destroyed  in  the 
area  A. 
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Weapons  Aimed  at  Small  Targets 

If  weapons  are  individually  aimed  at  a  target  whose  dimensions  are 
small  when  compared  to  the  errors  in  aiming,  i:he  variation  over  the  tar¬ 
get  in  the  probability  of  hitting  an  element  of  area  dxdy  can  be  neglected 
and  the  lethal  area  is  an  adequate  measure  of  destructiveness. 


The  bombing  errors  x  and  y  along  the  range  and  deflection  directions 
are  usually  assumed  to  be  normally  distributed  with  standard  deviations 

and  ay ,  respectively.  Then,  the  probability  of  hitting  a  target  element 

dxdy  is 

2  2 

Pt(x,  y)  dxdy  =  (2ff  a  a  )  1  EXP  {-  |  [—■  +  ^-j]}  dxdy  .  (II-9) 

y  x  y 


Near  the  target,  i.e.,  at  x  =  y  =  0, 

Pt(0,  0)  dxdy  -  ^  dxdy.  (11-10) 

Consequently,  the  probability  of  destroying  the  target  with  a  single 
weapon  is 

=  jj  Pt(x,  y)  f  U,  y)  dxdy 


_  |* [>  f(x,y)  dxdy  , 

JJ  2flO  a 
x  y 


(11-11) 


where  f(x,  y)  is  the  probability  density  function  of  hitting  the  target. 


Hence , 

p^1)  =  rf  f(x,  y)  dxdy  _  L 

k  <•’*)  2ir o  a  2w  a 

x  y  x  y 


(11-12) 


If  N  weapons  are  fired  independently  at  the  target,  the  probability  of 
killing  the  target  is 


(1 


L 


2770  o 
x  y 


(11-13) 


Since  by  assumption  L  is  small  compared  to  o  and  o  ,  equation  (11-13)  be- 

x  y 

comes 


(N)  ,  -NL/ lv, a  a 

p  1  -  e  x  y 


(11-14) 


Equation  (11-14)  is  derived  following  the  same  reasoning  used  to  obtain 
equation  (U-8). 
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Weapons  Aimed  at  Large  Targets 

When  the  target  dimensions  are  such  that  the  assumptions  for  small 
targets  are  no  longer  valid,  the  variation  over  the  target  of  the  prob¬ 
ability  of  hitting  an  area  element  dxdy  must  be  considered,  and  the  lethal 
area  concept  is  no  longer  useful. 

Again,  let  f(x,  y)  dxdy  be  the  probability  of  hitting  an  area  element, 
th'ir.,  the  probability  that  the  target  is  destroyed  by  a  single  weapon  is 


pk1}  =  JT  pt(x>  ■>  dxdy 


(11-15) 


and  further  progress  hinges  on  the  ability  to  evaluate  this  integral.  Once 
elevated*  ho  ever,  one  has,  once  more,  the  result  that  the  expected  number 


of  shots  required  to  destroy  the  target  is  i/P^\  and  the  probability  of 
destroying  the  target  with  N  weapons  is 


-  (1  -  P 


(1KN 
k  ’ 


(11-16) 


which 


if  P, 


(1) 


is  small,  can  be  written,  as  before, 


(11-17) 


When  Pt(x,  y)  is  constant  over  the  target  area,  equation  (11-15)  repre¬ 
sents  the  probability  of  hitting  the  target,  multiplied  by  p^,  and,  if  the 
probability  density  function  f(x,  y)  is  not  too  complex,  the  integration 
can  sometimes  be  carried  out.  This,  however,  occurs  seldom  in  practice, 
and  it  is  precisely  the  evaluation,  in  closed  form,  of  the  integral  type 
given  by  equation  (11-15)  that  has  been  the  primary  concern  of  the  re¬ 
searchers  involved  in  Task  301. 


Fragment  and  Blast  Sensitive  Targets 

A  fragment  and  blast  sensitive  target  is  one  in  which  the  major  dam¬ 
age  mechanism  is  due  to  fragmentation  and  blast  effects  rather  than  to  a 
direct  impact  by  the  weapon.  In  this  case  the  damage  function  p  (x,  y) 
is  the  kill  probability  for  a  warhead  given  that  it  has  detonated  at  (x,y); 
f(x,  y)  is  the  probability  density  function  for  the  warhead  detonating  at 
(x,  y). 

In  this  case  the  integration  is  carried  over  the  entire  xy  plane.  It 
is  with  integrations  of  this  type  that  the  investigation  reported  here  was 
primarily  concerned. 
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SECTION  III 


ANALYTICAL  EVALUATION  OF  AVERAGE  KILL  PROBABILITY 


The  evaluation  of  the  average  kill  probability  integral  given  by 
equation  (11-15)  has  formerly  been  accomplished  by  numerical  techniques; 
programs  written  by  Martin  Marietta  Company  are  available  for  this  purpose. 
Such  techniques  involve  random  samplings  of  the  damage  function  employing 
Monte  Carlo  procedures. 


In  order  to  investigate  the  nature  of  the  damage  functions,  regression 
analysis  was  applied  to  numerically  defined  damage  functions  generated  by 
computer  programs  developed  for  the  Air  Force  Armament  Laboratory  by  Martin 
Marietta  Company.  The  regressed  functions  had  to  possess  the  property  of 
having  high  coefficients  of  correlation  as  well  as  being  integrable  in 
closed  form  when  convoluted  with  the  probability  density  function  f(x,  y) 
which  was  assumed  bivariate  normal  with  zero  coefficient  of  correlation, 
i  .  e  .  , 


f(x,  y)  = 


2flo  o 
x  y 


1  X'^x 

EXP  f-  j  U~ 


y-Hy 

)  +  — -1 


2 

■)  ]  } 


(III-l) 


where  p  and  p  are  weapon  biases  due  to  aiming  and  ballistic  errors  in 
x  y 

the  range  and  deflection  directions,  respectively,  p  and  p  are  the  co¬ 
ordinates  of  the  weapon's  mean  point  of  impact  for  a  target  Centered  at 
the  origin  of  the  coordinate  system. 


The  damage  function  pt  (x,  y)  depends  upon  the  target  definition  as 
well  as  upon  the  following  weapon  parameters: 

1.  Height  of  Burst 

2.  Terminal  Velocity 

3.  Elevation  Attack  Angle 

4.  Fragment  Mass  or  Masses 

5.  Impact  Pattern  Dimensions 

The  numerically  defined  damage  functions  are  best  described  in  polar 
coodinates;  for  this  reason  it  was  necessary  to  transform  all  mathematical 
functions  to  this  system,  thereby  giving 


f  (  r  ,  8) 


1 


2ru  o 
x  y 


EXP 


{  -  \  t(I 


cos  9  - 


j.  \  ‘ 
X  ) 


+ 


r  sin  9  -  u  ~ 
( - - - -1) 


y 


]  } 


( III-2) 
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where  x  =  r  cos  0, 
and  y  =  r  sin  0. 


(III-3) 


2-rr  <° 

=  J  J  Pc(r.  9)  f<r.  0)  r  dr  d0 

o  o 

2ff  [a  (9)r2  +  0(0)r] 

=  T1  J  J  e  r  dr  d0  ,  (III-6) 

o  o 


where 

«<e>  =a(0)  -  + 

x  y 

P(0)  =  B(0)  +  u  cos  0  +  u  sin©  ,  (III-7) 

x  y 


and  T1  =  - -  EXP  [C-ip,3--^8!  . 

1  2ito  a  2  tc  2  y 

x  y 

The  complexity  of  the  functions  required  that  a  Maclaurin  series  expansion 
of  pt  (r,  9)  f  (r,  9)  be  generated  in  order  to  integrate  them  in  closed 

form.  This  was  accomplished  as  described  in  Appendix  I,  with  the  integra¬ 
tion  along  the  coordinate  axes  performed  from  0  to  rQ  and  from  0  to  9^; 

where  r  is  a  limiting  value  of  the  radius  where  the  regressed  equation 


7 


for  p.  (r,  8)  applies,  and  8  is  one  full  cycle  (360°  angle).  Both  r  and 

K  O  O 

8o  must  be  smaller  than  unity  to  speed  up  convergence.  The  following 
result  was  obtained: 

8o  ro 

=  J  J*  Pt(r,  9)  f  (r,  0)  rdrdB  = 
o  o 


~C  ■  2  '  2  8  r  2 


l-(j-2)  j+1 


.  i-1  r 

1  r  „  o 


(6  ) 


2  it  o  a 
x  y 


f— °.  °_  +  y  _ i.  f  y  f1')  — - 2 - 

L  2  +.^  i!  L.^0  V  [i-(j-2)](j+l) 


g°(i-j)r,je]}. 


(III-8) 


where 


,m-l. 


gmr,n0  ^  V  (m-l)r ,  (n-k)8Pr\k8+(  1  )g° (m-2)r(n-k)8 
k=0 


o 

p2r,ke}][l-sm(0)]+sin(0)sn(0),  for  m>0,  n>0,  (III-9) 

o 

Pmr,  n0  =  [2b  s  (2)  +  b.s  (1)  +  b.s  (0)]  sm  (1) 

A  n  z  n  J  in 

+  [4a  s  (2)  +  2a  s  (1)  +  (2a.  -  - oT^sCO)  ]s  (2) 

In  2n  3  za  z<j  n  ■*  m 

x  y 

+  fuxc°S  (^")  +  UyCOS  [  0^)1*]]  sn  (1) 

♦  ,2"  <2T--ir->  ( !*> 1  *. (2>'  (1II‘10> 

y  x 

for  m  =  1 ,  2;  n  >  0, 
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s\(j)  =  1  for  k  =  j 


s,  (\i)  =  0  for  k  /  j 


(III-ll) 


Typical  output  from  the  regression  program  is  given  in  Appendix  II. 
Future  Work\ 


The  Maclkirin  expansion  solution  is  not  satisfactory  due  to  its  slow 
convergence  to\P^.  Nevertheless,  the  average  kill  probabilities  evaluated 

following  this  procedure  are  more  reliable  than  the  numerically  computed 
ones  because  they  are  obtained  from  the  consideration  of  a  larger  number 
of  sample  points  and  also  because  of  the  upgrading  of  information  about 
the  nature  of  the  damage  functions  achieved  by  the  method  developed. 

To  overcome  the  ^disadvantages  of  slow  convergence,  the  following  pro¬ 
cedure  is  proposed  for\  future,  continued  research  in  this  field: 

I.  Evaluation  of  Damage  Functions 

R.  Snow  and  M.  Hyan  (Snow,  1968,  p.5)  has  assumed  the  damage 
function  p  (x,y)  of  equation  (3)  to  be  of  Gausian  form: 

2  2 

(Equation  A)  p  (x,  y\  =  D(x,  y)  =  D  exp  (-D  [^r-  +  ]) 

\  Rd)  R(2) 

Based  upon  this  assumptio^  they  have  developed  many  useful  closed 
form  relationships  for  calculating  fractional  coverage  to  areal 
targets. 

Preliminary  work  on  this  task\  indicates  that  the  symetry  suggested 
by  this  form  may  not  be  supported  by  the  numerical  functions,  par¬ 
ticularly  for  low  elevation  angles  of  attack.  Also,  for  detonations 
of  high  height  of  burst  (HOB)  which  are  out  of  the  blast  range  of 
the  target,  the  damage  function  takes  on  lower  values  near  the  origin 
of  coordinates  due  to  the  relatively  small  fragmentation  densities 
in  the  polar  zone  nearest  the  nose\of  most  ballistic  weapons. 

In  order  to  evaluate  the  validity  of\  the  Gausian  damage  function^ the 
following  study  is  proposed.  \ 

1)  Test  the  hypothesis  that  the  damage \function  is  of  the  form 
given  in  equation  A  against  the  alternative  that  it  is  not. 

2)  Study  modifications  of  the  form  given  in  equation  A  to  explain 
the  unsymmetrical  cases  described  above  and  test  the  hypothesis 
that  the  damage  function  is  of  the  form  given  by  equation  A 
against  the  alternative  that  it  is  of  the  modified  form. 
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II.  Optimization 


The  explicit  formulation  of  damage  functions  for  single  weapon 
single  target  element  cases  makes  possible  the  study  of  more 
complex  cases  of  multiple  weapon  attacks  on  multiple  element 
targets,  by  integrating  the  same  damage  function  with  different 
density  functions. 

One  problem  that  has  not  formerly  been  treated  is  that  of  opti¬ 
mizing  the  pattern  parameters  of  a  multiweapon  attack  on  a  dis¬ 
cretely  defined  target  complex.  It  is  proposed  that  optimization 
techniques  for  optimizating  these  pattern  parameters,  using 
damage  functions  described  in  Part  I  above,  be  studied  by  this 
task. 

Optimization  Algorithm 

During  the  summer  of  1969,  the  principal  investigator  collaborated 
with  other  L.S.U.  personnel  in  teaching  a  course  in  Systems  Optimization 
Theory  at  Eglin  Air  Force  Base.  In  part  as  a  result  of  this  experience, 
the  principal  investigator  developed  an  acceleration  algorithm  to  minimize 
a  convex  objective  function  subject  to  linear  constraints.  A  paper  de¬ 
scribing  the  algorithm  has  been  submitted  for  publication  in  the  Journal 
of  the  Operations  Research  Society  of  America  and  is  included  in  this  re¬ 
port  as  Appendix  III.  It  is  anticipated  that  the  algorithm  will  prove 
to  be  relevant  to  the  work  for  Task  301,  because  the  researchers  expect 
that  the  average  kill  probability  function  will  be  convex  and  that  the 
ballistic  constraints  will  be  expressible  as  linear  inequalities. 

Reference 


Snow,  R. ,  and  Ryan,  M. ,  A  Simplified  Weapons  Evaluation  Model  . 
Report  No.  RM-5677-PR,  The  Rand  Corporation,  December,  1968. 
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APPENDIX  I 

MATEHMATICAL  DEVELOPMENT 
Consider  tne  average  kill  probability  integral, 

OO 

=  If  ■f(x,ii)c#.7LJ!u 

.00  ^ 
zrTo* 

~  Jj  (r,e)  f(r,e)r<Hr&e  j 


O  o 


ex/3 


(^yi]  j 


rt ('x, 3)  =*■  11  (r,  e)  =  ex/”  / [A(a)r *  t&(e)rtC]}i 


where 


<e)  -  a.,©*  -ha. ^6  +  Q-3  ; 
/3fe)  -  b,ez  4-bx  G  +  i 

C,  ~  Constant: 


and 


ra  =  %*■  y *•  X=  r  C8*»  ©y  y  -  r -o**-  O  - 

exe 


Therefore, P ( ^  ©)  “  ZTTrr^a^  w  t  u  r  /  '  °  $ 

but  ^0-Hx^  =  r^^-Z^rcae+l^x1, 


and  6"  ©  “  H  <3 


-  MO*  r  rl  xm^q  -  2.  My  r  awe  ©  +  2. 


Hence, 


^  f( 


r  cpog  -  H* \ 

2  ,  (rat*.  ©-My  U 

tf>  / 

+  <-  >s  / 

[«7 

{^r  x  coo  *■  ©  -  2  Mx  i 

*vLr 

*  /we 2  ©  -  Z  r  aix/  0  + 
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+  M  -AM-t 

£(e)  =  Bfe)  +  H  X  &  +  U  tj  Ain-  G  , 

=  C  “  z:  HtC.Z  "  ir  Wy*-  CONSTANT. 


Then, 


//  7%  (rte)  f  (r^)  r&r  <6 G  - 

.  -  v*  Cx7r L^C.o'u-*-  +&(e')r+jJ 

-  (zrr  (Tic,  <r$ )  /  /  eL  r^^cffe- 

°  9 

x-i  a.  r°°  r2  f/?  ('©)0 

=  (2.7T  (TX  ^ <j)  e/  J  J  G  rcQrcfie 


Then, 


tfy)~'er=^ 

zTT^ 

J r f  ?°t  (r,  &)  { (r,e)  r  Ar  A  G  - 

nr  «, 

*Y  ff  <t  1>‘e)ri~'fe)rJ  rrfrdf-e  . 


Consider  a  Taylor  series  expansion  of  a  function  of  two  variables 
about  point  (VQ  }  0O), 


3^,e)  =  J 


^  fro  .  e0) 


^  =  T(r-rJ  j;  +  (e-So)  ^  ] 

If  T~0  —  0o  a  obtain  the  Maclaurin  expansion 

or  3  O,  e)  as  follows: 
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1 


9 


(r,  e)  -  £ 


i- o 


where  the  operator 


4  =  Cr£*e-£] 


Now  consider  the  function 


[j,(e)r*+JS(e)rJ  _  (f  g) 


=  e 


Where  (V,  ©)  =  (<*-,  +  GLfl  0  +  CL  3  ~  Ceo  6  ~  £<r^-  r' 

+  Cb,exftae  *-b3  +  Hx  c«o6  +  Nij  oaa ©)  r 
=  £*  >-  *  +  j3  r 


Let  •  a.5  -  Z<T„*  J 


b4=  Hx  ,  = 


Then, 


§£=Tm  =  + 


<*JS 

^r* 

ay 

<?r* 


-  ^r)  =  *  * 


=  ^(jr)  -  ^ 


rr5  ~  *  O  for  -rr  ^  3 

n  t 
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■AimmIMIMs 


^  -  fza.te  +  <zz-  cl  +  2.  +a.5zc**e 

+  tz-b4  **a.q  +  k5  Coo  ©)r 

-  [zat9  +  Q-2. -  fa.  4.-  0_5)  2.  Coo  ©  >o>wc  ©Jr*- 
+  L  2 1,  ©  f  t*.  -  t>4  a**-  0  +-  t>5  Coo  ©J  r 

A)  - 

~  L*- a-,  ~  (<*-  +  - a.s)  z  C~' a^z©-+  C6X5*o)Jr 

■h  C.  Z  b ,  *  CsO  8  -  by  A*^  ®J  r 

=  L  Z  a-i  -  z.  (o*a  -  o^s)  Cl-  Z  r* 

•f  £  2  2> ,  ~  k  4  Coo  ©  *  A1*1  ©Jr 

*  5  ^o>  *  4^e  c*06>)]  r* 

+■  L  i>4.  aaA  ©  ~  4^5  c®°  ©Jr 
=  £  8  (a.*-  «-*)  W©  co©  ©  ]  r  * 

+  L  fc>  4  <ua  ©  -  ^5-  coo  ©Jr 

J^.  =  tf?4S)=  LBfa.,- 

■h  £.  fc>  4.  (L^o  9  f-  ^5-  ©*■*-  ©  J  r 

=  <Pc<-e)  =  c“®Jr‘ 

+  £.-  i>4  a*^  ©  +■  co^  ©-Jr* 

-  <^4,©)=  L‘3z(cx.+-  <ZS)G-Z  W*©)jr* 

j,  *1 

Hh  L-  O^o  ©  ~  r 


§5;  =  7e>=  L/28(aA-ac)(A^e<*»e')]rg' 

+  L  ^4  av>-  ©  -  ^5  c^°  ©J  r 

if??  =  <5^.=  L 118  (*A-d;)(i- 2.  «+«.*&) ]rz 
<Je  * 

+-  £  fc>  4.  Cj®o  ©  +  £>5-  a*a-  ©Jr 


§p=<n*e>  =  CeM-'V^+-06-z^eJr* 

+•  £  £>4  ©  r  £>^.  «m  &  3  t* 

for  ~n  =  4,  &,  >z,  IG,  ZO,  ... 

£>r  *>,.=  L-z‘~*(o.a-  *s)0-z~~*o)3r* 
ae"  '  _  _ 

-h  L  -  £>4  C©«  O  -  b  y.  AM.  &  Jf" 


n=  <o,  IO,  I  A~,  13,  22, 


©jsp 


=  tf^e)  m  L~z*  (ol4’*s)  cex>ejrz 

+  t"t>  4  0  +  c«o  ©Jr 

yi*  S,  9,  J3,  17,  Zl,  ... 


6^Ue)  =  ^©Ctfoejr1 

+  £  t> 4  q-w^.  e  -  bs  c*o  e]  r 
for  v-3,  7>  ">  2  3, 
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Further, 


"  ^e)  -  C_0  *  [l.  £**  Vo.*-  o.s)(\-Z*J*-6Y]tx 

+  Lb4  ©  +  by  a***.  ©3 

for  -yi  =■  A- 1  Q>,  Q ,  I  O  t  IZ,  1 4,  ... 

=  ^0  *  [Lzn  (a4- as)jiU& cocgJ  rz 

f-  £  b  4  /w>-  e  -  Ls  c^c  ej  r"] 

for  *=  <3^7,  ?  //,  /3,  ... 

'»  jp  <0,  o)  *  O 

0=^0,  o')  -  ^6  To)  -  Jb3tk  4 

CP(lLr)  Co, O)  -  2-  <*  fo)  *  £  C(X  3  +CL  4) 

^3r)  C°>0)  *  ^4r)(0js-  B  ^r)  CO)  *  O 

0(e)  C°,o)  =  tfr*®)  Co,o)  =  ...  -  ^e)  Co,0)=0 


Computing  Mixed  Partial  Derivatives  cfy3  Cr,e) 

a 


^~r~a~  =  ^r,e)  =  2  L  2-ol,g  +  at  -  Col4-(Xs)  z  too e *»*.&]  r 

+  L  Z  b,  ©  +  b^  -  b  4  W.  0  +  b5  <?J 


6>V 


r  ^  -  ^Vr,  ^  L  Z  a. ,  ~  Z.  Co.  4~  Q-f  )  C I  ~  Z 

+  L  Z  b,  -  b  4  9  -  b5  <Q>/c  ©] 


/a***  4  ©X]  r 
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a^V  .  s  1  Z,U*Ca4’As)  eao  sir 

^rae"  (r-*«  r  ,  T? 

+  £  b4  ^  e  -  b5  <^oej  > 

for  -yi  -  -3 ,  £ ,  7 ,  ^ ,  ll i  ■  ■ 

*'*>  -  -p  --  C-tf^  |2[£hJ’(M-il^|-z“‘,e)lr 

I  +■  L  t>4  C-®»  0  +*  bj-  .a>-A-  wj  J 

\  .  __4  (o ,  8,  IO  1 2.  ■  ■ 

?!f  ,  r£ir  e,  =  2  Cz ». ®  - “■»  *  u*~  ^ Zc^e ™ e3 

<9*r}e  '  ’ 

a  V  arf>  .  -  2£?«a>-2Ca+-aPO-^^i®0 


a _ d£L  - 

dr^e" 


_  =  c-o”*1  iiie(a+-zi)*~B<»*t>Pi 


for  -*«-  =  -3,  f,  7,  %  t>,  •  •  • 


^  -  c-1)^  \^Lt^  C0.4.-  dsYi-^^l} 

j^e1*  VUr.ney  > 

^  ^  *  s  Q.  t  r\  1  *> 


for  „J  4,  <..6,10,11,... 


?  <£  -  -O  =0  for  '  Vt'1 


Finally: 


,0  (0,0)  *  ba  +  W 

^r,e)  ’ 


CO,  o')  -  2  b,-  bi 
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fp  Co,o)  -(-!)  *  C-bs)  for  n- 3,5,7%  II, 


■n-4 


Jp  Co  o)  =  r-0  *  (b*)  £°r  n  =  4,  (a,  8,  I O,  iz,  ■■■ 

v  (r,no) 


f%rr.»,(0'O)  =  2a- 

=  Z[g.a,-2(a.  4-as)l  =  4(a.,-a.4  +  0-e) 


'Ur  no) 


(0,0)  =  O  for  3,5,7, 9,  //,... 


^2r, no)  (o,o)  =  (r I )  4  [2£2.fn  ,Va4-^)3ls^0  *  2n(°-4-<£) 

for  /l=  •4-,  & ,  © ,  ^O/  ••• 


^r,n»)  ~  ^  for  ^ 


2  3,  all.yx  . 


Summary 

_ ^  fo,o)  -  o 

— ^<>>(0,0)  -  b3  +  k* 

— ^  r°>0)  a  2ro.3^a4) 

— ^  Z^c)  M)  "  O  for  --7  *  3 


^  =  °  for  «  e| 
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*%•)  (°>0)  =  b*  +  b* 

iPi.r.zs')  ~  2  b, 

tr(r  ne>)  (o,o)  =  f-.)V  (b^  for  n  =  3,S.  7.%  II,  ... 

^<r,Ae)  55  (“*)  *  C  b4)  for  r»  =  4. <o, 8,10,1*,.-. 

(°>°)  =  2^-2 

<P(ar,ae) (o.o)  =  4(fl.,-a4  +  as> 4a.  ♦£»)  2Y«-4-^ 

(0,0^=  O  for  n  *  3,ff,7fS,  II,... 
f(,r,fl9)  (0,0>  =  (-»Yr  2n  (*-♦-  a5)  for  vt=4,C,,8,/0,»2,... 
=  °  for  »2?-  n>'* 


Let  n«)  (0,0")  =  J^„r>  rt® 


The  formulas  for  ^»r,n«  can  all  be  condensed  into  the  following 
expression: 


'Cr.nB  =  [  t  2t,An(«.)+bzArt0)  +y3*n(o)]a'^U) 

f  [Ad,  A*n(-d  +  ZOl^a yv(/)^C2.a3+  0.4.  f  Aj)An  *,(£■) 

+  C{t>4COQ^^^r)  +  by  CsoLC  'sl')tt]  1  A/m  0) 

+  frra4-a,)W^ff)l^^]  tt  £7-^(0)] 

j=3  ^ 


Where  >0^  =/ 

^  6)  =  O 


*=J 

*  /j 


DISCRETE  UNIT 
STEP  FUNCTION 
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V*' 


Sr-  &fft  *3 ft 


Sir'- S' ft  *3  ft' 

frr'&r  ft  +  Sr  ftr  +  Sr  ftr  *  ffSr 
=  S„ ft  *  Zfrftr  +§ftr 

S-,/S,rft  'Skrftr  *  ZStrftr  +  ZSr^r'Srftr  + 
=  Sir  ft  +3f,rftr  +3S-rftr  *fftr 
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r 


f  i 


Sir  foX  3&r  t„  +  3  %,  fr 

+  3  for  for*  3S-<  for  +  Sr  for  *  $fs. 

~-  for  *r*  *%,  f,r  +  Gfor  P3r  +*fo  for  *  ftsr 


% 


KS  7Ui,n  *  ^  fo~.*Xr  fi*J  % -»vF,r 


tz) 


*-  +(~^fo  *L.,)r  +  ff-~r 


for 


But  /**  *  O  for  -m.  ^3.  Therefore, 

9  1*9  r 


o  _ 


O  /  >*~H  * 


t~r  "  fo^rfo  *■  (  r'hL,)r  Kr  £“ 


Similarly, 


fo;  +  *  ^*ht*-*e  foe 

*■■■*  G'-A  ^  S' foe 


But 


n 


_  o  for 


we 


>  / 


i 


The 


refore,^ 


»  O  for 


■w© 


>t 


i  / 
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■a  — 


To  obtain  mixed  partial  derivatives,  consider, 

h*m  *•***''■• 

fr.ro  =  %•*'  *  f.fr.o  +  %  fro  +  f  fr.ro 
=  fi.fr  *  *  %  Tr.o  *  ffr.ro 


fr.ro  ‘  fo.fr  +  fcofr.o  * 


+  9*ff.*e  *  30 

“f..n 

+  •*&.»$•  +  3 

$r,  *® 

+  frfr,3* 

no  ~  *3-rofr 

+  (  ‘)  §■<•*-}»  fro 

,  *(D 

C^Cn-3L)0  y^X.® 

+  .  .  . 

+  %  7? 

*($ 

^  r, 

for  Tt  Jl 

Therefore , 


SC) 

K*0 


^  (-«-&)  e  y° r.  If  e 


for  -yt  i  / 


‘tf-nT*  hrJ 

%r  hrofr  *  frfr..1  *  l%f~*  fKrJ 
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+  §r.e  fr,e  +  Hr,oTr,m  +  frrfr.zo^ 

*  %e^zr.z  '  %'o'fte.e  *  3  $lr.  zfi^ 

~  Llr.t.fr*  *9r.eK,ei'  f'Tr.zel 
+  £,}«.nr  +  21.5ir,e  +  fKr.J 

ftr,ie  "*  ^'^r.(A-rteJr.8+^*^5r.(«-J)»,r««  ' 

+  (~,)$r,a1°r.<+-«»*($%Tf.«»] 

+  ^ne  ^er  +  ^  ^n.,)e  lar.o  +  ^  fxr.ze*  '  '  ' 


The-efore, 


for  't'Z  "2-  I 


Finally, 


7*r,*ie 


for  V1  ^  I 


In  general. 


o  of  vn-i\  0  o 

^wt-'V  +  ^  '  '  ^(w»-a)r  ^a  f 


for  «?l 
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From  which 


~  ^  ^r,  Is 

^  ^5,1. ,W  U-L\e  ffrkJ 


+  t  '  rt(*-h)9  ti*r,9to 

for 


Summary 


0) 


(  vrt  *  71  =  O') 


&)  £r  =  &-)r  «■*  + 


$6*-z)r  'tf’xr 


for 


vrr  -  / 


<a>£.  =°  £or 


^  =  JE  (*)  { Cw-.v,  *»  -  *>«  *  • 

<7(->*r.  w©  «,  _  v 

”  °  /  ^  _  l\  o  o  ) 

+  l  i  /  S  fw-*)r,  Cti  -  Me  'Gzr,  he  J 

for _ tYI  ^  I  ,  tT  ^  I 

Where 

-.  j[2b,^n(z^^b3tQ,rt0Vb3^(o)]>^/m6) 

+  [4a,A  a)f  2atAn(i) 

+  ^as'£c7  '  zcr^^n.^^  ^ 
+  &©*  (li  +  1  Am  6) 

+  $£*(53^  ~‘z<rx*~)t**  jl^  (a)l 

IT  £i-aw0>} 

•>=j  J 
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Now,  the  Maclaurin  expansion  of  ^  C^t&)  is  given  by, 

c.&,e)  *  £  °2 

c/  i,Q  i  l 


But 


Cr &*e£li9(o,o).gr  r‘*-6‘)  jfrJr,  r*l 

+-  r'“‘e'  *  ••• 

re>1"  +  0 -)fite 


-  ^  /*'\  o  ^  ^t-j5  i 

bzo  4  ®  * 


Consequently, 


Now  , 


( r,  e)  rdtrcQe 


can  be  computed. 


For  practical  and  obvious  reasons,  the  integration  along  the  T coordinate 
will  be  performed  from  O  to  t"0  ,  where  is  some  limiting  value  radius 
where  the  regressed  equation  for  ter,*)  applies.  Also,  the  integration 
on  6  will  be  carried  out  from  O to  ©„  where  ©„  is  one  full  cycle  (360° 
angle).  Both  and  0,,  must  be  smaller  than  unity,  to  speed  up  convergence. 
This  is  accomplished  by  proper  scaling  of  the  coordinate  system. 


Therefore 


*  r. 

f  I $(r’  e)  r^r<0e  - 
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&*  r. 


i-G-i> 


i 


e  <&r<Qe 


1=0 


<‘-jKj 


e 


*-<!-«) 

r. 


<•-6'-*) 


.? 


©  l  (Q  & 


Observe  that,  consistent  with  notation. 
Therefore, 

<9.  r0 

]  f Cfe)r  <9r  <Qe  - 


">t»  ~  1  ~J  and  72  =7 


o»  ;  j+.  *1 

a  <  jL  [s  (l)  - - <=*-• 

*'  fr-^Ci+O  T‘<-i K.J®J 

-  1 

u.i( 

fciKj.] 


<-  S  ,t 

is  i 


i-  Ci- *)  i*  • 


i\  r.1  ’  -  e». 


Ci-ti-vld+ri 


Ci-j)r ,  j® 
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and 


But 


9  °  =  O  =$?  for  -yrz  =  o  Ci-j) 
&  n.  B 


Then, 


'"o 


// 


•  o 


rc&re#e  — 


--  ^  +i  .t 


t-i 


i  1  r. 


<>i 


jso 


i-fj-l)  j^i 

ea 

<r<v-jV, 


J 


Finally: 


©o  r0 


Z/^T 


«)/ 6)  a)r&f<$Q  - 


o  o 


[C-1M**-  i*?} 

> 

z,'tr<r%<r^ 


&0  n 


+ 


£ 

im 


-<i-z)  i+l 


b-(j-zj]  6*o 


Where 


o 


1  )  *«}]  +J>m  (°) 


for 


^  O,  'yl  ^  O 
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C \n9-  t,An(0+  b2/>n(t)  +  b3  AnCo)J  Am  (l) 

+  (0  +Za3LxnLt)  +  tza3- £rt  ~  )  An (o\)  (z) 

+  ^Mx  Cd* it)  +  Wij c**i(  t')ttJ  J  ^ y*i  f () 

+  <**(?”)]  ft) 

for  -wt  —  1,2.)  -rt  -O 
AJ»  (j)  =  /  for  J,  =j 
^  =  O  for  Jr  *,/ 


and 


APPENDIX  II 

RESULTS  OF  REGRESSION  ANALYSIS 


Curves  of  the  form  given  in  equations  21  and  22  were  fit  to  numerica1- 
ly  defined  damage  functions  for  weapon-target  case  152. 

The  damage  functions  used  are  listed  in  Table  II-I  along  with  the 
corresponding  goodness-of-fit  ratios  (GOF)  and  sums  of  the  squares  of 
the  deviations  of  the  curves  from  the  observed  data  (E(y-y)^) •  The  good- 
ness-of-fit  ratio  is  defined  as  the  ratio  of  the  variance  explained  by  the 
explicit  relationship  to  the  variance  of  the  numerically  defined  function 
(the  number  one  represents  perfect  correlation). 

For  high  height  of  burst  (HOB)  the  values  from  the  numerical  functions 
are  very  small  and  almost  randomly  scattered.  Curves  fit  to  such  data  will 
have  small  E(y-y)  since  the  deviations  from  the  zero  base  are  small  to  be¬ 
gin  with;  the  same  curves  may  have  very  poor  goodness-of-fit  since  they  ex¬ 
plain  very  little  of  the  variance  present  in  the  data.  It  is  possible 
that  there  is  no  relationship  between  such  data  points  (complete  randomness) ; 
small  values  of  2(y-y)“  are  still  possible. 

The  ninety-degree  elevation  angle  of  attack  (ANG)  cases  can  be  fit  with 
symmetrical  surfaces.  The  goodness-of-fit  ratios  for  these  cases  are  larger 
than  those  for  lower  elevation  angles  of  attack. 
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TABLE  II-I.  SUMMARY  OF  RESULTS  OF  REGRESSION  ANALYSIS 
(Weapon-Target  Case  152) 


Terminal  Conditions 


Regression  Curves 


HOB 

ANG 

VEL 

GOF 

L(y-y)2 

1  .40 

45.00 

700.00 

0.933629 

0.39671E  02 

1 .40 

45.00 

1000.00 

0.874630 

0.75235E  02 

1.40 

45.00 

1300.00 

0.955912 

0.26579E  02 

1.90 

70.00 

700.00 

0.958644 

0.24414E  02 

1.90 

70.00 

1000.00 

0.962145 

0.22360E  02 

1.90 

70.00 

1300.00 

0.960273 

0.23510E  02 

2.10 

90.00 

700.00 

0.997513 

0.15603E  01 

2.10 

90.00 

1000.00 

0.995336 

0.26659E  01 

2.10 

90.00 

1300.00 

0.994256 

0.31494E  01 

20.00 

45.00 

700.00 

0.795667 

0.13674E  02 

20.00 

45.00 

1000.00 

0.786459 

0.16873E  02 

20.00 

45.00 

1300.00 

0.777634 

0.20675E  02 

20.00 

70.00 

700.00 

0.767283 

0.15696E  02 

20.00 

70.00 

1000.00 

0., 821736 

0.16674E  02 

20.00 

70.00 

1300.00 

0.826441 

0.21226E  02 

20.00 

90.00 

700.00 

0.807875 

0.29011E  02 

20.00 

90.00 

1000.00 

0.837392 

0.27062E  02 

20.00 

90.00 

1300.00 

0.892173 

0.23002E  02 

40.00 

45.00 

700.00 

0.663889 

0.14383E  01 

40.00 

45.00 

1000.00 

0.653946 

0.18524E  01 

40.00 

45.00 

1300.00 

0.650797 

0.23304E  01 

40.00 

70.00 

700.00 

0.618090 

0.17297E  01 

40.00 

70.00 

1000.00 

0.669085 

0.21498E  01 

40.00 

70.00 

1300.00 

0.709352 

0.26624E  01 

40.00 

90.00 

700.00 

0.972755 

0.20575E  01 

40.00 

90.00 

1000.00 

0.967214 

0.28898E  01 

40.00 

90.00 

1300.00 

0.961912 

0.38102E  01 

80.00 

45.00 

700.00 

0.571082 

0.75003E  01 

80.00 

45.00 

1000.00 

0.539971 

0.10604E  00 

80.00 

45.00 

1300.00 

0.272525 

0.21975E  00 

80.00 

70.00 

700.00 

0.513242 

0.10193E  00 
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TABLE  II -I  (CONCLUDED) 


HOB 

ANG 

VEL 

GOF 

E(y-y)2 

80.00 

70.00 

1000.00 

0.544083 

0.13480E  00 

80.00 

70.00 

1300.00 

0.604326 

0.17633E  00 

80.00 

90.00 

700.00 

0.945677 

0.10548E  01 

80.00 

90.00 

1000.00 

0.950072 

0.12186E  01 

80.00 

90.00 

1300.00 

0.956756 

0.13884E  01 

120.00 

45.00 

700.00 

0.514104 

0.11593E  01 

120.00 

45.00 

1000.00 

0.435045 

0.18436E  01 

120.00 

45.00 

1300.00 

0.429930 

0.23478E  01 

120.00 

70.00 

700.00 

0.521283 

0.14438E  01 

120.00 

70.00 

700.00 

0.521283 

0.14438E  01 

120.00 

70.00 

1000.00 

0.529390 

0.18597E  01 

120.00 

70.00 

1300.00 

0.549753 

0.24629E  01 

120.00 

90.00 

700.00 

0.924258 

0.56311E  00 

120.00 

90.00 

1000.00 

0.928711 

0.68378E  00 

120.00 

90.00 

1300.00 

0.932487 

0.81306E  00 

These  results  and  the  actual  regression  coefficients  are  available 
from  the  researchers  in  punch  card  coded  form. 
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APPENDIX  III 

ACCELERATION  ALGORITHM  TO  MINIMIZE  A  CONVEX  OBJECTIVE 
FUNCTION  SUBJECT  TO  LINEAR  CONSTRAINTS 
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ABSTRACT 


An  acceleration  algorithm  to  minimize  a  convex  objective  function 
(one  that  is  never  underestimated  by  a  linear  interpolation  between  two 
points)  subject  to  linear  constraints  is  presented  in  detail  and  an  ex¬ 
ample  problem  is  given.  It  reduces  the  number  of  changes  in  the  state  set 
and  the  number  of  iterations  required  by  the  general  differential  algorithm 
for  linear  constraints  in  converging  to  the  minimum.  It  can  also  be  used 
to  maximize  a  concave  objective  function  (one  whose  negative  is  convex)  sub¬ 
ject  to  linear  constraints. 


34 
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INTRODUCTION 


An  objective  function  y  (x)  Is  convex  if  it  is  never  underestimated 
by  a  linear  interpolation  between  two  points,  x^  and  x2 .  That  is,  for 
every  a  satisfying 

o  <  a  1  (1) 

it  is  true  that 

y  +  [l-<v]  x2)  *  nry^)  +  [l-<v]  y(x2).  (2) 

The  adjective  "strictly"  is  added  if  the  function  is  never  exactly  equal 
to  the  linear  interpolation  between  the  two  points,  i.e.,  if  the  <  sign 
is  replaced  by  the  <  sign.  A  concave  function  is  one  whose  negative  is 
convex  (Zukhovitskiy  and  Avdeyeva,  Reference  2). 

If  a  non-linear  function  is  to  be  minimized  subject  to  linear  con¬ 
straints,  i.e.,  a  problem  of  the  form: 


Min.  y  (x) ,  where  x  =  I 

*N 


(3) 


subject  to. 


x  so  for  n  ■  1 ,  . . . ,  N 


and  N  (4) 

Z.  \n  Xn  *  V  k  =  X*  •••’  K’ 

n=l 

The  following  differential  algorithm,  given  by  Wilde  and  Beightler 
can  be  used  (Reference  1,  pp.  65-66): 

1.  Let  v,  be  the  most  negative  decision  derivative  and  v.  the  most 

i  n 

positive  decision  derivative  for  which  the  corresponding 
decision  variable  d^  is  positive. 

2.  If  there  are  no  negative  v  ,  set  v  to  zero.  If  all  positive  v 

mi  n 

have  corresponding  d  equal  to  zero,  set  v,  to  zero. 

n  h 

3.  If  both  v^  =  0  and  v^  =  0,  a  stationary  point  has  been  found. 

Notice  that  when  both  v,  and  v,  are  zero,  the  Kuhn-Tucker 

i  n 

necessary  conditions,  i.e.,  non-negativity  and  complementary 
slackness,  are  satisfied. 
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4. 


5. 


6. 


Tf  either  v.  or  v  is  not 

i  h 


zero , 


compute 


V  =  +  vh  (5) 

If  V  <  0,  increase  d^.  If  V  >  0,  decrease  d^ ,  holding  all  other 
decision  variables  constant,  but  permitting  the  state  variables 
to  readjust,  until 

I.  Some  state  variable,  say  s  becomes  zero,  or 

P 

II.  becomes  zero;  where  r  =  i  if  V  X  0,  and  r  =  h  if  V  >  0, 

or 


III.  d.  becomes  zero, 
n 

7.  In  case  (I),  d^  replaces  s^  in  the  state  set;  in  cases  (II)  and 
(III)  no  change  in  the  set  of  states  occurs. 

8.  Clearly,  for  any  change,  dd^,  the  value  of  the  objective  function 
y  always  decreases  since 


Ay 


3d  <  0. 
r 


(6) 


The  greatest  difficulty  in  the  application  of  the  preceding  algorithm 
occurs  in  step  7,  case  (I),  when  the  variables  d^_  and  s^  are  interchanged, 
i.e.,  when  a  change  in  the  state  set  becomes  necessary,  because  a  new  set 
of  constrained  derivatives  must  be  computed  as  follows: 


=  -.sy.  = 

*d 

r 


^2- 

ad. 


k=l 


kr 


_dy_ 

ds, 


;  r  =  1,...,N  (7) 


where,  v  = 
r 


Sd 


is  the  constrained  partial  derivative  of  the  objective 


function  y  with  respect  to  the  decision  variable  d  ; 

ay 

~  ~~  is  the  unconstrained  partial  derivative  of  the  objective  function 
3d 

r 

y  with  respect  to  the  decision  variable  d^_; 

“V  is  the  unconstrained  partial  derivative  of  the  objective  function 
k 

y  with  respect  to  the  state  variable  s^;  and 
is  a  constant  coefficient  obtained  from  the  following  linear 

equation  relating  the  state  variable  to  the  decision  variables 


d  ,  n=l , 
n 


,  N, 

=  8, 


N 

-  y 
n=  1 


o',  d  ,  k  =  1,  —  ,  K; 
kn  n 


(8) 


6^  is  a  constant  term  in  the  expression  for  s^ . 
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Acceleration  Algorithm 

When  the  objective  function  is  convex,  the  following  minimization 
algorithm  can  be  used  to  reduce  the  number  of  changes  in  the  state  set 
and  to  decrease  the  number  of  iterations  required  for  convergence.  The 
same  algorithm  can  also  be  used  to  maximize  concave  objectives  subject  to 
linear  constraints. 

I.  Compute  minimum  point  x°  for  unconstrained  objective  function. 

II .  Check  the  constraint  set  to  ascertain: 

1. 


III. 


IV. 


V. 


2. 


1. 


2. 


If  all  the_constraints  are  loose  at  x° ,  the  minimizing 
policy  is  x*  =  x° .  STOP. 

If  one  or  more  constraints  are  tight  or  violated  at  x° , 
CONTINUE. 

Let  the  slack  variables  of  the  constraints  that  are  tight 
and/or  violated  at  x°  be  in  the  decision  set  d. 

Let  the  slack  variables  for  the  constraints  that  are  loose 
at  x°  be  in  the  state  set  s. 

3.  Complete  the  state  and  decision  sets  with  the  structural 
variables  of  the  problem. 

Set  the  decision  variables  to  zero, 

d  =  5.  (9) 

Use  the  differential  algorithm  for  a  non-linear  objective  function 
subject  to  linear  constraints  to  find  the  minimizing  policy  and 
verify  that  the  point  obtained  satisfies  the  sufficiency  con¬ 
ditions  for  a  minimum. 


The  algorithm  guarantees  that  no  changes  in  the  state  set  will 
occur  provided  that  all  the  loose  constraints  at  x°  remain  loose 
at  x.*  Since  the  decision  variables  are  set  to  zero  initially, 
they  can  only  increase  in  value  to  satisfy  non-negativity.  That 
means  that  all  the  slack  variables  for  the  tight  and  violated 
constraints  at  x°  will  be  zero  at  the  beginning  of  the  algorithm, 
thus  making  these  constraints  tight  initially.  Since  the  decisions 
are  manipulated  at  will  and  cannot  go  negative,  the  initially 
tight  constraints  can  either  remain  tight,  if  the  corresponding 
slacks  remain  at  zero,  or  go  loose,  if  the  corresponding  slacks 
increase,  and  unless  one  or  more  of  the  loose  constraints  at  x° 
is  tight  at  x  *,  no  changes  in  the  state  set  will  be  required  and 
the  number  of  iterations  for  convergence  to  the  minimum  is  con¬ 
siderably  reduced. 


Illustrative  Problem  and  Conclusions 


Application  of  the  acceleration  algorithm  will  be  illustrated  with 
the  following  problem  involving  a  quadratic  (convex)  objective  function: 

Min.  y  =  2x^  -  2x^2  +  2 -  6x^  (10) 

subject  to  - 


xl’  x2  —  0; 

(11) 

3x^+4x2  ±  6; 

(12) 

■x^+4x2  <_  2. 

(13) 

Step  I.  To  compute  the  minimum  point  for  the  unconstrained  objective 
function  set. 


ay 

—  =  4x,  -  2x0  -6=0, 
ox,  1  2 


3y 

3x„ 


=  4x2  -  2x^  =  0  , 


(14) 


from  which 


(15) 


The  differential  quadratic  form  is  positive  definite,  that  is. 


. -T 

JX 


H 


ax 


(3x1> 


=  4  (aXi 


+  3 


(16) 


for  dx^,  3x2  f  0. 

This  guarantees  that  x0  is  in  fact  a  global  minimum. 
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_ _ 


Step  II.  A  check  of  the  constraint  set  reveals  that  constraint 
(12)  is  violated  and  constraint  (13)  is  tight  at  x° . 

Step  III.  Therefore,  after  the  introduction  of  slack  variables,  the 
problem  becomes: 


2  2 
Min.  y  =  23^  -  2s1s2  +  2s2  - 

subject  to  s^,  s2>  d^,  d2  >_  0; 


3s2  +  4s2  +  =  6; 


and  -s^  +  4s2  +  d2  =  2' 


with  d. 


X- j  d«  ~  xa >  ®i  “  xi *  —  xo • 


Solving  for  s.^  and  s2  in  terms  of  d^  and  d2  obtain 


81  *  1  "  *  4  dl  “  4 
32  =  4  “  (  16  dl  +  16  V ' 


Because  of  the  form  of  equation  (8) , 


81  “  1;  “ll  =  4 5  a12 


1 

4’ 


B2  =  4;  a21  “  16  ;  a22  =  16. 


Therefore,  from  equation  (7), 


v=«JL.  =  iX__(0l  +  a  &-) 

V1  6d1  3d1  '  11  3s1  °21  3s2  ’ 


and  v 


-  it  (52  +  13  dl  -  17  d2); 

J§2_  -  15L_  _  (a  12-  +  a  iZ_  N 

2  “  6d2  3d2  (a12  3Sl  “22  3s2  ) 

-65  H>8  ~  17  dl+  37  d2>- 


Step  IV.  Set  dx  -  d2 


0. 


(17) 

(18) 

(19) 

(20) 
(21) 

(22) 


(23) 

(24) 
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Step  V.  (1-1) 


Having  ■  d2  “  0| 
from  equation  (20)  , 

from  equation  (21) , 

from  equation  (23), 

from  equation  (24)  , 


obtain. 


-  l; 


S2  *  4  ’ 


3 

4 ' 
13 

'l  *  16 


17 
64 - 


(25) 


(2-1) 

(3-1) 

(4-1) 

(5-1) 

(6-1) 


The  differential  algorithm  requires, 

17  n 

V  “V  =»  —  —  -  •  v  =  11 

i  2  64  ’  h 

Both  v.  and  v,  are  not  zero,  therefore 
l  h 

the  minimum  has  not  yet  been  found. 

c-  17 

Since  v|  =  v2  =  comPute 

V  =  V.  +  V.  “  -  77. 
i  h  64 


Because  V  <  0,  increase  di  =  d2 . 

The  increase  in  d2  is  constrained  by: 


Case  (1) 


s^  +  Ask2  —  0  f°r  ^ 


1,...,  K, 


but 

iSk2 
Therefore , 

,  i  min 

Ad,, 


Ad2  =  -  ak2  Ad2 , 


fe) 

3re’  /  3  \ 

m±n  i  !i<_  1  L-\.  4 

\2  >0  '  ak2  '  \  ji/ 


Case  (II) 


6v2 

V2  +  7T 


Ad„ 


0. 


Therefore , 


Ad 


2 


(26) 


(27) 


(28) 

(29) 

(30) 

(31) 

(32) 
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Case  (III)  Since  d^  is  initially  set  at  zero  and  is  to  be 
increased,  this  case  does  not  apply. 


From  the  results  above,  Ad^  =  ^ 


(33) 


(7-1)  No  state  change  was  needed,  and  the  new  values  of  the 
decision  and  state  variables,  and  of  the  constrained 
derivatives  are: 


dl  =  ° 


remains  constant; 


,  ,  ,  68  _  68 

d2  °  +  37  “  37  ’ 

_4  .  68  54 

31  “  1  “  16  '  37  37 ’ 

8  _  1  ,_3v  68  .  15. 

2  4  '16''  37  37’ 


17  X  M  =  12 
64  '37  37  ’ 


v  .  _  1Z+  (  37  68  „  o 

2  16  ^  64  '  37  - 

198 

The  value  of  y  at  this  point  is  y  =  -  . 


(34) 


(35) 


The  results  can  be  displayed  in  a  tableau  with  the  following  format 


dl 

d2 

S3 1 

all 

“12 

8i 

“21 

a22 

S2 

— 

V1 

V2 

y 
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Therefore, 


0 

SI 

4 

4 

54 

16 

16 

37 

1 

3 

15 

16 

16 

37 

12 

-  198 

0 

37 

37 

Now  the  differential  algorithm  is  applied  all  over  again. 

(1-2)  The  new  values  of  decision  and  state  variables  as  well  as  of  the 
decision  derivatives  are  given  in  the  tableau. 

(2-2)  v.  =  0  and  v  =  0  (36) 

l  h 

(3-2)  Since  both  v.  =  0  and  v^  =  0,  the  minimum  has  been  found  at 


1 

SI  = 

54 

;*=  4  = 

37 

15 

(37) 

37 

and  y* 

(x*)  = 

198 

37  • 

(38) 

The  problem  converged  to  the  solution  in  one  iteration  and  no  changes 
in  the  state  set  were  necessary. 

The  same  problem  was  solved  by  Wilde  and  Beightler  (Reference  1,  pp. 
76-78)  using  the  differential  algorithm  alone,  without  acceleration.  Three 
iterations  involving  two  changes  in  the  state  set  were  required  to  con¬ 
verge  to  the  minimum. 
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